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摘要 : GNSS 信号 在 卫星 发 射 以 及 接收 机 接收 的 传输 过 程 中 存在 设备 时 延 。 设 备 时 延 与 钟 差 、 
模糊 度 等 密切 相关 ， 很 难 将 其 与 它们 完全 分 开 。 通 常 ， 根 据 实际 需要 将 设备 时 延 与 其 相关 参数 组 
合 形成 不 同 的 相对 设备 时 延 。 如 果 能 更 好 地 将 这 些 相 对 设备 时 延 进行 分 类 和 区 分 ， 对 高 精度 的 
GNSS 应 用 具有 重要 的 参考 意义 。 分 别 介 绍 了 目前 基于 设备 时 延 提 出 的 几 个 常用 偏差 : 码 间 仿 
差 、 未 校准 的 相位 时 延 、 频 间 钟 差 偏差 和 系统 间 偏 差 ， 详 细 说 明了 它们 的 概念 、 产 生 原 因 、 解 算 
策略 及 相关 研究 进展 。 
关 键 词 : 码 间 偏差 ， 未 校准 的 相位 时 延 ， 频 间 钟 差 偏差 ， 系 统 间 偏差 
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1 引 言 


在 全 球 卫 星 导 航 系 统 (global navigation satellite system, GNSS) 的 发 展 历程 中 ， 美 国 
和 俄罗斯 的 卫星 导航 系统 不 断 地 发 展 并 进行 现代 化 改善 ， 其 他 国家 和 地 区 同时 也 和 争 相 建设 
自己 的 卫星 系统 并 发 展 相 关 产 业 ， 不 同 卫星 导航 系统 之 间 相 互 融 合 与 借鉴 ， 极 大 地 推动 了 
GNSS 技术 的 进步 。 因 此 ，GNSS 的 服务 能 力 不 断 提升 ， 达 到 了 质 的 飞跃 ; 4 GNSS 技术 相 
关 的 应 用 领域 不 断 扩 大 ， 已 成 为 人 们 生活 中 不 可 缺少 的 一 部 分 。 在 GNSS 的 应 用 日 益 广泛 
的 背景 下 ， 研 究 者 和 用 户 越 来 越 不 满足 于 现在 的 成 绩 ， 追 求 更 高 质量 的 数据 、 更 好 的 数据 处 
里 方 法， 以 及 获取 更 精密 的 产品 为 其 主要 目标 之 一 。 在 导航 定位 相关 方面 的 应 用 中 ， 能 够 更 
好 地 改正 、 估 算 甚 至 完全 消除 GNSS 信和 号 传输 过 程 中 的 各 种 误差 ， 是 获取 高 精度 产品 的 先 
决 条 件 。 其 中 ，GNSS 设备 引起 的 时 间 延 人 运 (简称 为 设备 时 延 ) 是 主要 误差 源 之 一 。 

GNSS 的 设备 时 延 (或 称 之 为 硬件 延迟 ) 是 指 GNSS 信和 号 在 卫星 和 接收 机 传输 过 程 中 受 


到 设备 通道 环境 (如 其 中 的 数字 滤波 器 ) 影响 产生 的 时 间 延 迟 。 其 中 ， 将 信号 从 产生 点 到 达 
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卫星 天 线 相位 中 心 ， 经 过 设备 通道 的 时 间 延 迟 称 为 卫星 发 射 设备 时 延 ， 从 接收 机 天 线 相位 中 
心 到 真正 信号 处 理 点 的 时 间 延 迟 称 为 接收 机 设备 时 间 延 迟 。GNSS 伪 距 和 载波 相位 观测 方程 
中 ， 设 备 时 延 通常 作为 方程 中 的 附加 项 ， 包 括 接 收 机 和 卫星 相关 的 设备 时 延 参数 。 而 在 实际 
应 用 中 ， 卫 星 和 接收 机 的 设备 时 延 与 相应 的 钟 差 高 度 相 关 ， 相 位 设备 时 延 也 与 模糊 度 有 关 ， 
不 能 完全 分 离 ， 同 时 卫星 与 接收 机 的 设备 时 延 也 无 法 完全 区 分 开 。 因 此 ， 一 般 在 数据 处 理 过 
程 中 ， 设 备 时 延 作为 相对 偏差 进行 解 算 。 表 1 列 出 目前 常 采用 的 处 理 与 设备 时 延 有 关 偏 差 的 


4 种 方法 "。 
表 1 设备 时 延 相关 偏差 的 处 理 方法 
方法 方法 介绍 
消去 法 将 偏差 从 观测 方程 中 完全 移 除 ， 一 般 使 用 差分 的 方法 


参数 法 将 偏差 作为 观测 方程 中 的 参数 ， 同 其 他 参数 一 起 估计 处 理 
广播 修正 法 ”使 用 其 他 方法 获得 偏差 并 实时 播发 给 用 户 ， 该 时 延 在 一 定时 间 段 内 稳定 
校准 法 事先 根据 其 他 来 源 进行 偏差 的 解 算 并 作为 常数 使 用 ， 该 时 延 一 般 具 有 长 期 稳定 性 


在 不 同 的 应 用 中 ， 需 要 考虑 的 设备 时 延 略 有 差异 ， 根 据 实 际 需 要 ， 对 设备 时 延 有 关 的 
偏差 进行 分 类 研究 。 在 单个 GNSS 系统 中 ， 涉 及 到 两 个 及 以 上 信号 的 伪 距 处 理 时 ， 由 于 不 
同 信号 /频率 在 通道 中 的 传输 时 间 不 一 样 ， 需 要 考虑 信号 间 相 对 伪 距 设备 时 延 ， 即 码 间 偏差 
(differential code bias, DCB)。 当 进行 模糊 度 解 算 时 ， 为 了 固定 整 周 模糊 度 ， 出 现 了 关于 未 
校准 的 相位 时 延 (uncalibrated phase delays, UPD) 的 研究 。 当 使 用 多 个 频率 的 载波 信号 时 ， 
不 同 频率 组 合 下 的 卫星 钟 差 存在 一 定 差异 ， 需 要 考虑 频 间 钟 差 偏差 (inter-frequency clock 
bias, IFCB)。 在 进行 多 个 系统 联合 处 理 时 ， 就 产生 了 系统 偏差 (inter-system bias, ISB) 这 一 
概念 。 

随 着 国际 上 对 设备 时 延 相 关 偏 差 的 研究 日 益 增 多 、 应 用 趋 于 广泛 ， 国 际 GNSS 服务 
组 织 (The International GNSS Service, IGS) 成 立 了 偏差 和 校准 工作 组 (The IGS Bias and 
Calibration Working Group, BCWG)， 该 工作 组 主要 是 对 GNSS 偏差 领域 进行 相关 的 研究 。 
它 主 要 目的 是 探讨 目前 存在 的 设备 时 延 定义 与 影响 ， 给 出 适当 、 一 致 的 处 理 规则 ， 为 用 户 提 
供 相 应 的 高 质量 产品 。 同 时 ， 关 注 不 断 出 现 的 新 GNSS 系统 和 新 信号 ， 不 断 更 新 设备 时 延 
的 定义 、 增 加 新 内 容 ， 完 善 提 供 的 产品 服务 。 

本 文 将 对 与 设备 时 延 相 关 的 GNSS 偏差 进行 详细 介绍 。 


2 FSB) mize 


由 于 测 距 码 在 传输 或 者 接收 过 程 中 产生 的 设备 时 延 不 一 致 ， 在 相同 时 刻 不 同 频率 间或 者 
同 频 不 同 测 距 码 间 出 现 的 设备 时 延 偏差 ， 即 为 DCB”, DCB 是 两 个 信号 间 的 设备 时 延 之 
差 ， 值 的 大 小 受信 号 通道 的 影响 ， 可 达 几 纳 秒 甚至 十 几 纳 秒 ”"。 根 据 信号 机 制 ， 可 将 DOB 
分 为 频 内 DCB 和 频 间 DCB， 前 者 是 指 在 同一 个 载波 上 因 测 距 码 不 同 造 成 的 设备 时 延 偏 
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差 ， 而 后 者 指 由 于 载波 频率 不 同 产生 的 设备 时 延 偏差 ”。 
接收 机 DCB 和 卫星 DOB”. 

DCB 是 反 演 电离 层 总 电子 含量 (total electron content, TEC) 的 主要 误差 源 ”"。 早 在 
1991 Œ, Coco 等 人 在 对 GPS 系统 的 电离 层 延 迟 精 度 进行 确定 时 ， 发 现在 工 波段 之 间 存 在 
设备 时 延 相关 的 偏差 ， 并 认为 在 确定 电离 层 总 电子 量 时 必须 将 该 偏差 从 GPS 观测 中 去 除 ”。 
在 GNSS 的 精密 单 点 定位 、 时 间 同 步 等 方面 的 应 用 中 ， 也 需要 考虑 DCB”。 所 以 当 涉 及 到 
两 个 及 以 上 信号 的 伪 距 观测 的 处 理 时 ， 需 要 考虑 到 DCB"”。 此 外 ， 由 于 接收 机 DCB 与 接 
收 机 的 设备 性 能 密切 相关 ， 对 接收 机 DCB 数值 进行 长 期 分 析 有 利于 研究 接收 机 的 变化 ”。 
随 着 高 精度 GNSS 的 广泛 应 用 ，DCB INVER AA BRR, DCB 的 研究 也 趋 于 完善 。 关 于 
DCB 的 处 理 ， 频 内 DCB AME DCB 有 所 区 别 。 

关于 频 内 DCB 的 处 理 ， 可 通过 伪 距 观测 值 直 接 做 差 解 算 ”。 因 为 两 个 码 信号 调和 
同 的 载波 频率 上 ， 可 以 忽略 电离 层 的 影响 ， 伪 距 观 测 方程 见 式 (1): 


于 延迟 产生 的 位 置 不 同 ， 又 分 为 


在 相 


= 


P=p+c(dt, — dé) +(d,-@)+T+I+ep, . (1) 


将 两 个 码 P,, Po 的 伪 距 观测 值 直接 做 差 可 得 DCB”， 如 公式 (2): 


P, — Pp = (dri — dy2) + (d — di) + ep. , (2) 


其 中 ， 尸 表示 伪 距 观测 值 ，p 表示 为 接收 机 和 卫星 间 的 几何 距离 ，dt., dts 分 别 表 示 接 收 机 
和 卫星 钟 差 ， 必 , ds 分 别 表 示 接 收 机 、 卫 星 的 伪 距 硬件 延迟 ，7 为 对 流 层 延迟 ， 了 为 电离 层 
延迟 ，epl, eps 为 伪 距 噪声 改正 。 由 于 伪 距 观测 值 形成 的 观测 噪声 较 大 ， 无 法 完全 忽略 。 可 
将 上 式 观测 值 在 一 天 内 取 平 均值 ， 削 弱 噪 声 ， 得 到 综合 DCB 观测 值 。 

频 内 DCB 参数 在 一 段 时 间 内 (如 一 个 月 ) 比较 稳定 ， 可 将 其 看 作 一 个 常数 ”。 欧 洲 
定 轨 中 心 (Center for Orbit Determination in Europe, CODE) 2010 年 开始 用 上 述 方法 解 
算 GPS 和 GLONASS 系统 DOB 产品 并 提供 给 用 户 。 图 1 为 该 中 心 2019 年 1 月 卫星 频 
内 DCB 和 接收 机 频 内 DCB 的 月 均值 ， 图 2 展示 了 2019 年 1 月 G02, Gos KÆ ALRT, 
COCO 测 站 的 频 内 DCB 的 时 间 序 列 图 。 

频 间 DOB 参数 可 采用 事先 标定 和 软件 估算 两 种 方式 得 到 ”"。 卫 星 和 接收 机 在 出 厂 前 ， 
rr 通常 会 标定 设备 时 延 ， 用 户 可 直接 使 用 ” 。 随 着 设备 使 用 时 间 的 增加 ， 卫 星 和 接收 机 老化 
~ 导致 硬件 性 能 会 发 生变 化 ， 同 时 在 外 界 环境 等 多 种 因素 的 影响 下 ， 实 际 值 与 标定 值 产 生 差 
异 ”。 Coco 等 学 者 将 预 发射 校 准 值 与 估计 值 进行 比较 ， 结 果 表 明 四 对 卫星 中 有 两 对 显示 出 
极 好 的 一 致 性 ， 而 另外 两 对 差异 显著 ， 认 为 需要 谨慎 对 待 事先 标定 值 ”"。 因 此 更 多 的 研究 者 
和 用 户 利用 GNSS 的 实测 数据 对 频 间 DCB 进行 解 算 ， 并 且 这 些 解 算 的 数据 可 对 标定 的 频 
间 DCB 进行 监测 ”。 
目前 ， 关 于 频 间 DCB 的 估算 方法 有 两 种 : (1) 在 电离 层 TEC 建 模 时 ， 将 频 间 DCB ES 
数 与 其 他 参数 进行 计算 ”™ ，(2) 采用 经 验 或 已 有 的 电 高 层 模 型 进行 电离 层 延迟 改正 ， 然 后 
再 计算 频 间 DCB 参数 “。 方 法 二 中 解 算 的 DOB 值 的 质量 很 大 程度 上 受到 选取 的 电离 层 
模型 精度 的 影响 一 。 但 当选 定 模型 后 ， 可 以 对 任意 系统 的 DCB 进行 计算 ， 并 且 可 随时 解 
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VE: 这 4 幅 图 从 左 到 右 依次 为 GPS/GLONASS 的 P1-C1 #4, P2-C2 码 间 卫星 频 内 DOB 和 接收 机 频 内 
DC 的 月 均值 示意 图 ，* 代表 GPS 系统 ，o 代表 GLONASS 系统 。 图 中 每 个 点 代表 该 系统 下 每 个 卫星 或 者 
测 站 对 应 的 频 内 DCB. 


1 2019 年 1 月 GPS/GLONASS 的 频 内 DCB™ 


1 = 
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图 2”2019 年 1 月 G02, Gos 卫星 和 ALRT, COCO Wàt P1-C1 DCB 的 时 间 序 列 中 


算 DCB， 所 以 方法 二 较为 方便 快捷 ， 极 大 地 提高 了 DCB 解 算 效率 ”。 德 国 宇航 中 心 采用 
了 方法 二 ， 即 利用 全 球 电 离 层 图 (global ionospheric map, GIM) 直接 扣除 电离 层 TEC， 获 
得 卫星 和 接收 机 的 综合 DCB， 然 后 基于 零 均值 基准 约束 方法 分 离 卫 星 和 接收 机 DCB 参 
数 一 。 采 用 该 方法 的 不 足 之 处 是 在 电离 层 变 化 显著 的 地 区 (如 赤道 )，DCB 的 精度 会 随 之 下 
降 ”。 而 方法 一 的 DCB 是 与 电离 层 模型 系数 同时 解 算 ， 结 果 会 受到 电离 层 解 算 精 度 的 影 
响 ， 所 以 需要 选取 分 布 均匀 并 且 数 量 合理 的 GNSS 测 站 一 。 
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此 外 ， 目 前 在 进行 电离 层 建 模 时 ， 假 设 接收 机 DCB 日 变化 稳定 ， 数 据 处 理 时 一 天 计算 
一 个 参数 空 。 对 于 采用 码 分 多 址 技术 的 卫星 系统 来 说 ， 如 GPS (global positioning system), 
Galileo、 北 斗 等 ， 这 样 的 处 理 方 式 是 合理 的 。 但 是 GLONASS (global navigation satellite 
system) 卫星 系统 与 前 面 三 者 不 同 ， 它 采用 的 是 频 分 多 址 技术 来 分 离 不 同 卫 星 传输 信号 。 
在 这 种 复 用 技术 中 ， 为 不 同 卫 星 的 载波 分 配 了 频带 L1 和 L2 内 的 相 令 频率。 不同 的 载波 频率 
言 号 在 接收 通道 中 处 理 的 方式 不 同 ， 延 迟 也 会 因 频 率 的 不 同 有 差异 ”; 所 以 , 将 GLONASS 
系统 中 因 载 波 频率 不 同 而 在 传输 或 者 接收 通道 中 产生 的 设备 时 延 的 差 值 称 之 为 频 间 偏差 
(inter-frequency bias, IFB)"". GLONASS 卫星 的 IFB 的 值 相差 极 大 ， 处 理 方法 上 与 频 间 
DCB 略 有 不 同 ， 需 要 将 每 个 接收 机 通道 中 的 IFB 解 算出 来 。 目前 IGS 分 析 中 心 提 供 的 
DCB 中 不 包含 IFB。 


3 未 校准 的 相位 时 延 


es 与 伪 距 类 似 ， 载 波 信号 的 发 射 和 接收 也 会 受到 各 种 偏差 的 影响 ”。 在 初始 相位 、 卫 星 
ES 和 接收 机 设备 时 延 及 其 他 因素 的 影响 下 ， 非 差 处 理 算出 的 载波 相位 模糊 度 包含 了 其 他 偏差 ， 
pt 被 模糊 度 吸收 的 这 部 分 相位 偏差 表示 为 UPD”. HP RABBI RET NH, UPD 参数 
= 的 整数 部 分 与 模糊 度 完 全 耦合 ， 无 法 分 离 ， 一 般 只 能 测定 不 足 载波 信号 一 周 的 部 分 ， 将 模糊 
To) 度 不 为 整 周 的 小 数 部 分 称 为 小 数 周 偏 差 (fractional cycle biases, FCB)”。 一 些 学 者 有 时 会 
O 直接 用 UPD 代替 FOB 进行 表述 。 
Cy 精密 单 点 定位 (precise point positioning, PPP) 作为 非 差 数据 处 理 中 常见 的 定位 方式 ， 
© 在 相关 应 用 中 常常 要 考虑 UPD。PPP 主要 是 利用 精密 轨道 和 精密 钟 差 等 产品 ， 在 综合 考虑 
N 各 项 误差 改正 的 基础 上 ， 采 用 合理 的 参数 估计 方法 (最 小 二 乘法 或 卡尔 曼 滤 波 法 等 )， 利 用 
= 单 台 GNSS 接收 机 实现 绝对 定位 的 技术 ”。 由 于 PPP 采用 的 是 非 差 观测 值 ， 许 多 误差 不 能 
通过 组 成 差分 方程 来 减弱 或 者 消除 ， 因 此 各 种 误差 源 都 必须 加 以 考虑 。 载 波 相位 模糊 度 的 
固定 ， 可 以 提高 PPP 的 精度 和 收敛 速度 呈 。PPP 技术 在 获取 大 气 延 迟 、 动 态 精密 定位 、 地 
壳 运 动 及 海洋 潮 位 监测 、 精 密 授时 、 精 准 农业 领域 等 方面 均 得 到 了 很 好 的 应 用 ”。 而 更 好 、 
更 快 地 固定 模糊 度 是 推动 PPP 发 展 的 重要 因素 。 


但 是 直接 计算 出 来 的 模糊 度 不 具有 整数 特性 ， 是 一 个 浮 点 解 ， 即 存在 UPD， 主 要 产生 
原因 有 三 个 ”: (1) 由 于 卫星 和 接收 机 产生 的 载波 信号 并 非 从 零 相位 开始 ， 并 且 初始 相位 未 
知 ， 无 法 将 其 从 模糊 度 中 分 离 。(2) 不 同 的 信号 在 卫星 端 和 接收 机 端 传输 中 ， 由 于 信号 的 差 
异 ， 会 产生 不 同 信号 间 设 备 时 延 ， 而 该 延迟 会 与 模糊 度 耦 合 。(3) PPP 使 用 的 精密 钟 差 产品 
主要 是 IGS 或 者 是 与 IGS 解 算 方式 相近 的 精密 产品 ， 它 们 在 计算 钟 差 时 一 般 使 用 无 电离 层 
组 合 形式 并 以 伪 距 观测 方程 的 钟 差 为 基准 ;由 于 其 精度 不 高 ， 只 有 厘米 级 ， 在 载波 相位 观测 
方程 中 ， 使 用 该 钟 差 基准 时 ， 模 糊 度 参数 也 会 吸收 部 分 钟 差 。 

与 模糊 度 浮 点 解 的 PPP 相 比 ， 基 于 模糊 度 整 数 解 的 PPP 优势 主要 有 两 个 :提高 定位 的 
精度 和 缩短 收敛 时 间 。 李 星星 等 人 ”指出 浮 点 解 的 PPP 要 达到 厘米 甚至 毫米 级 精度 ， 往 往 
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需要 数 小 时 以 上 的 观测 时 间 ， 并 且 其 定位 精度 和 可 靠 性 不 及 双 差 固定 解 。 近 年 来 ， 国 际 上 的 
研究 重点 转向 PPP 固定 解 ， 其 关键 是 恢复 非 差 模糊 度 的 整数 特性 ， 因 此 提出 了 分 离 FCB 
的 方法 。 
越 来 越 多 的 研究 开始 用 计算 FCB 的 方法 进行 模糊 度 固定 ， 其 中 最 常见 的 方法 有 星 间 单 
差 估计 法 和 整数 钟 法 。 星 间 单 差 估计 法 主要 由 Cabor jth, ZJ Ge 等 人 ™ 根据 他 提出 的 
方法 做 出 部 分 改进 : 首先 是 利用 多 个 (100 左右 ) IGS 地 面 观测 网 算出 非 差 模糊 度 的 浮 点 解 ， 
然后 在 卫星 间 做 差 消 除 接收 机 端的 FCB; 再 利用 MW 组 合 估计 宽 埠 模糊 度 ， 将 其 取 整 得 
到 宽 埠 FCB 改正 数 ， 再 将 宽 埠 模糊 度 带 入 无 电离 层 组 合 中 ， 固 定 窗 埠 模糊 度 ， 分 离 窗 埠 
FCB 改正 数 。 而 整数 钟 法 则 是 计算 非 差 宽 埠 的 PCB AUER, FIA HAE FCB KERK E 
定 宽 埠 模糊 度 ， 然 后 估计 包含 了 窗 埠 FOB 的 卫星 钟 差 ”。 当 获得 FCB 后 ， 直 接 对 无 电离 
层 组 合 下 的 宽 埠 和 窄 巷 模糊 度 进 行 固 定 ， 得 到 PPP 的 模糊 度 整数 解 ， 获 得 高 精度 的 结果 ， 
可 更 好 地 应 用 到 实际 中 。 
目前 武汉 大 学 PRIDE (Positioning Racers to Image and Decipher the Earth) 课题 组 提 
HET GPS C1W/C2W/LIC/L2C 的 卫星 FCB， 该 数据 是 每 颗 星 每 个 信号 一 天 一 个 值 ”。 
图 3 为 2020 年 1 月 1 日 所 有 GPS 卫星 4 种 信号 的 FCB 值 。 由 此 可 见 ， 在 同一 天 中 ， 不 同 
TÆ Af SH FCB 差异 明显 。 图 4 表示 为 G01 和 G10 两 颗 卫 星 在 2020 年 1 月 份 的 时 
间 序 列 图 ， 可 见 在 该 时 间 段 内 同一 卫星 同一 信号 变化 较 稳 定 。 法 国 CNES 分 析 中 心 也 提供 
了 宽 埠 TCB， 主 要 是 利用 整数 钟 法 解 算得 到 。 


G01 G10 G20 G30 


GPS PRN 


3 2020481 A 1 AFE GPS PÆ C1W/C2W/L1C/L2C 的 FCB™ 
与 只 有 GPS-FCB 方法 相 比 ， 利 用 多 个 系统 对 不 同系 统 的 FOB 值 估计 并 进行 应 用 ， 能 
够 获得 更 快 更 好 的 定位 结果 。 武 汉 大 学 测绘 学 院 也 在 进行 多 GNSS 系统 FCB 产品 的 研究 ， 
包括 GPS, Galileo, BDS 和 QZSS 等 系统 号。 目前 ， 关 于 FCB 的 研究 虽然 趋 于 成 熟 ， 但 还 
在 发 展 阶段 ， 提 供 的 产品 也 在 逐步 完善 中 。 
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4 2020 年 1 月 G01, G10 PÆ C1W/C2W/L1C/L2C FCB 的 时 间 序 列 四 


4 ” 频 间 钟 差 偏差 


美国 的 GPS. BRINE Galileo、 日 本 的 QZSS (quasi-zenith satellite system) 和 中 国 的 北 
斗 均 提供 多 频 信号 服务 ， 多 频 信号 的 加 入 进一步 推动 了 GNSS 的 发 展 。 使 用 多 频 信号 提高 
了 定位 精度 和 模糊 度 的 解 算 速度 。 相 比 于 双 频 信号 ， 三 频 信 号 具有 增加 可 观测 的 波长 、 降 


低 噪声 、 减 少 电离 层 影响 等 优势 ”。 利 月 


三 频 信 号 可 进行 周 跳 探测 和 修复 、 模 糊 度 固定 等 


方面 的 应 用 ， 推 动 了 高 精度 导航 定位 的 发 展 。 但 受 了 卫星、 接收 机 硬件 延迟 、 空 间 环境 等 影 


啊 ， 采 用 不 同 频率 观测 、 观 测 


组 合 解 算 的 卫星 钟 差 存在 差异 ” 。Montenbrucket 等 学 者 在 研 


J GPS Block IF 卫星 时 ， 发 现 不 同 频率 无 电离 层 组 合 估计 得 到 的 钟 差 之 间 存 在 一 定 的 差 


异 ， 并 将 其 定义 为 频 间 钟 差 人 


ia 22 (inter-frequency clock bias, IFCB), [FCB 最 大 可 以 达到 


15 ecm”。 并 且 他 们 还 认为 IFCB 的 变化 周期 与 卫星 运行 时 相对 于 太阳 高 度 角 的 变化 相关 
E, Li 等 学 者 在 其 基础 上 对 PRN25 和 PRN01 卫星 的 IFCB 进行 计算 ， 并 提出 了 解 算 频 


间 相位 偏差 的 方法 ， 用 于 估算 两 个 Block UF 卫星 的 IFCB， 研 究 发 现 IFOB 的 变化 对 不 
同 频率 上 无 电离 层 组 合 的 PPP 收敛 时 间 与 定位 精度 均 产 生 影响 所 。 
在 三 频 处 理 中 ，GNSS 三 频 观 测 可 以 组 合成 两 组 无 电离 层 组 合 的 形式 为 : 


Ermis =p + (cdt, 


Liris =p + (cdt, 


dy, ») 
dr, ») 


(cdt? 
(cdt? 


di2) +T + epi 


àir2 N12 +T + er 


di2) + (bria — drma = b12 F d*;9)4 


Prr, =p + (cdt, 


dy, 5) 


Lr, 二 0 十 (cdt, + drs) 


A15 N1,5 + T + EL2 


(cdt? 


dis) +T ep 


(cdt 


dis) + (brs — dry, — 01,5 + d°1,5)4 


其 中 1, 2, 5 RRR, Pre, Lip 分 别 表示 无 电离 层 组 合 下 的 伪 距 和 载波 相位 观测 值 ，p, dr, 
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与 公式 (2) 中 相同 ，b., O 分 别 为 接收 机 、 卫 星 的 相位 硬件 延迟 ，》 为 无 电离 组 合 载波 相 
位 的 波长 ，N 为 相位 模糊 度 ，er 为 相位 噪声 改正 。 当 采用 两 组 无 电离 层 组 合 进行 卫星 钟 差 
解 算 时 ， 参 数 化 的 卫星 和 接收 机 钟 差 会 吸收 对 应 的 硬件 偏差 ， 所 以 会 产生 两 组 卫星 、 接 收 机 
钟 差 。 将 cdt + de, 。 为 频率 1、2 组 合 下 解 算 的 接收 机 钟 差 用 Fi RR cdt + di2 为 频 
率 1、2 组 合 下 解 算 的 卫星 钟 差 用 fs; 表示 ; cdh 十 di, ,为 频率 1、5 组 合 下 解 算 的 接收 机 
钟 差 用 571s 表示 ，cdts 十 dis 为 频率 1、5 组 合 下 解 算 的 卫星 钟 差 用 dis 表示 。 但 同时 解 
算 两 组 接收 机 和 卫星 产品 ， 不 仅 增加 了 时 间 成 本 也 加 剧 了 计算 负担 四 。 为 了 提高 解 算 效 率 ， 
提出 了 一 种 消除 钟 差 产品 不 一 致 性 的 方法 ， 即 求解 IFCB。 主 要 思路 是 解 算 其 中 一 个 频率 
组 合 的 钟 差 并 将 其 作为 基准 ， 其 他 的 钟 差 可 用 该 基准 加 上 FCB 表示 中， 即 ， 


we (4) 


IFCB, = Or1,5 = Or1,2 
假设 相位 设备 时 延 被 模糊 度 吸收 , 将 不 同 的 频率 组 合 (1/2 和 1/5) 做 差 ， 其 中 接收 机 


和 卫星 的 几何 距离 以 及 对 流 层 延迟 会 被 消去 ， 仅 留 下 模糊 度 项 和 IFCB， 而 盆 距 的 仅 留 下 
IFCB， 具 体 如 公式 (5) Atm: 


ea ekg 


: (5) 
Lprrizs = Ling — Lirias = IF CB, — IFCB +AipNig — MN,2N1,2 + Ers 


RP, Mo, Nis 分 别 为 频率 1/2 和 频率 1/5 组 合 下 吸收 了 相位 设备 时 延 的 模糊 度 。 研 究 
认为 接收 机 频 间 钟 偏差 为 一 常数 ， 卫 星 频 间 钟 偏差 除了 常数 部 分 之 外 还 有 对 应 的 变化 部 


nl, 


TFCBs =6°+ifchs ， (6) 


其 中 ，6s 表示 IFCBs 的 变化 部 分 ，ifcbs 为 IFCBs 的 常数 部 分 。 卫 星 频 间 钟 偏差 的 变化 部 
分 可 通过 相位 观测 方程 进行 解 算 。 当 历 元 间 没 有 周 跳 发 生 时 ， 将 公式 (5) 的 载波 相位 观测 方 
程 进行 历 元 间作 差 ， 整 周 模糊 度 被 消除 ， 随 之 消除 的 还 有 常数 FOB, Fl LFCBS 的 常数 部 
分 。 当 利用 单个 测 站 ， 多 个 历 元 (如 个 ) 进行 计算 时 ， 每 个 历 元 的 变化 部 分 为 : 


k 
al 
Òn = Lprrizs(m) 一 天 >， 了 DTFl25 ; (7) 
k=1 


被 模糊 度 吸 收 的 部 分 可 以 通过 伪 距 观测 方程 进行 解 算 。 但 直接 对 历 元 进行 求 和 取 平 均 ， 得 到 
的 是 卫星 和 接收 机 常数 部 分 的 和 。 为 了 得 到 ifcbs 部分， 消除 接收 机 部 分 ， 可 选择 一 个 参考 
卫星 进行 相应 的 求解 。 
卫星 IFCB 的 变化 部 分 上 共有 时 变性 ， 可 建 模 为 与 太阳 -卫星 -地 球 角 相关 的 周期 函数 ， 接 
WHLIFCB 可 以 看 作 一 个 常数 ”。 所 以 IFCB 可 以 提前 被 预 估 出 来 ， 应 用 于 相关 领域 。 
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5 系统 间 偏 差 


系统 间 偏 差 是 多 个 卫星 导 
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阮 系 统 进 行 联合 导航 、 定 位 等 应 用 时 ， 必 须要 考虑 的 偏 


关中 下， 因为 不 同 的 GNSS 系统 采用 的 坐标 和 时 间 基准 不 同 ， 不 同系 统 信号 的 结构 、 体 制 


等 有 差异 号 ， 当 不 同 GNSS 信 


号 在 多 模 接 收 机 通道 中 传输 时 ， 会 产生 接收 机 设备 相关 的 时 


延 偏差 ”。 所 以 利用 多 系统 GNSS 数据 进行 融合 解 算 时 ， 需 要 考虑 这 些 系统 性 的 偏差 ”。 
RUA ISB 是 一 个 卫星 系统 (如 北斗 ) 观测 值 与 参考 卫星 系统 (如 GPS) 的 观测 值 放 一 起 


处 理 时 ， 所 需 考虑 的 一 个 改正 ”。 在 实际 应 用 中 ， 由 于 系统 间 坐标 基准 的 差异 主要 体现 在 


卫星 位 置 上 ， 在 多 模 GNSS 


数据 处 理 时 一 般 转换 成 相同 的 坐标 系 ”， 可 以 忽略 坐标 基准 差 


对 ISB MRI. PO, ISB 主要 包括 了 时 间 基 准 之 差 和 接收 机 设备 时 延 ™ 


随 着 多 个 GNSS 系统 的 


系统 间 偏 差 研 究 相 应 产生 。 人 解决 ISB 问题， 可 让 多 GNSS 系统 更 好 地 联合 应 用 ， 可 用 卫星 


出 现 ， 很 多 学 者 对 多 模 GNSS 导航 定位 的 方法 进行 了 研究 空 一 ， 


数量 剧 增 ， 从 而 增加 了 同一 时 刻 同 一 地 区 观测 到 的 卫星 数 ” ， 卫 星 分 布 的 几何 结构 得 到 了 


优化 ， 降 低 了 位 置 精 度 因子 ， 


进行 导航 定位 


有 利于 在 环境 恶劣 地 区 (能 够 观测 到 的 卫星 数量 较 少 的 地 区 ) 


目前 ， 国 内 外 关于 ISB 的 研究 ， 主 要 包括 ISB 的 处 理 方法 研究 ，1SB 的 来 源 及 特性 
HH, ISB 建 模 预测 研究 等 。 


无 论 在 伪 距 还 是 在 载波 相位 观测 方程 中 ，75 均 会 存在 


但 在 不 同 的 情况 下 ，7SB 的 研究 与 定义 咯 有 差异 。 当 多 模 GNSS 联合 网 解 时 ， 整 网 内 只 估 


算 一 个 基准 钟 差 。 但 不 同 的 信号 频率 和 不 同 的 卫星 系统 会 产生 不 同 的 设备 时 延 ， 而 这 些 时 延 
会 被 接收 机 钟 差 吸收 ， 因 此 不 同 卫星 系统 的 接收 机 钟 差 实际 上 存在 差异 ” ， 所 以 需要 在 不 


同系 统 中 引入 一 个 偏差 参数 。 
时 时 ， 一 般 选取 GPS 的 接收 机 钟 差 作为 两 个 系统 的 接收 机 钟 差 ， 用 cd + dra 表示， 将 引 


入 的 偏差 参数 用 ISB 表示 ， 


例如 在 利用 GPS (用 G 表示 ) 和 BDS (用 C 表示 ) 进行 联合 处 


则 两 系统 联合 解 算 的 无 电离 层 组 合 观 测 方 程 为 


Prp = po 


( 
Lon = p° + (cdt; + dra 
( 


(cdt® + d°) +T + ep, 
(edt? + d°) + (bra — dra — bS + dF) HANS +T + ero 


) 
) 
PS, = pF + (cdt, d:c) 
也 六 = p? + (cdt, drc) 


(cdt + d°) + ISB +T + eps , (8) 
(cdt? + d°) + ISB + (bro — drc — bF + do) 


ANF +T + ers 


ISB 表示 为 系统 C 相对 于 系统 G 的 系统 间 偏 差 。 如 果 考 虑 不 同 


a 


星系 统 间 时 间 基 准 的 差 


FR, ISB 中 还 存在 一 个 常数 偏差 参数 ”， 用 DD 表示 ， 则 TSB 可 表示 为 : 


ER ISB 公式 表示 的 是 采 


+5 7 NA KM 


ISB=(do—da)tD . (9) 


1054) Z HELA GNSS AAW ABI mz. 4 GLONASS TE 
星系 统 ， 采 用 频 分 多 址 技术 来 区 分 不 同 的 卫星 ， 它 与 其 他 GNSS 系统 间 的 ISB 相 比 ， 还 包 
。 因 此 GLONASS 的 每 颗 卫 星 都 要 计算 出 1 个 15B 值 ， 或 者 将 
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每 颗 卫 星 的 频 间 偏 差 去 除 。 值 得 注意 的 是 ， 当 ISB 参数 与 其 他 轨道 产品 联合 解 算 时 会 产生 
秩 亏 ， 因 此 需要 引入 一 个 额外 的 约束 条 件 ” 。 目 前 有 两 种 约束 方法 ; 一 是 选择 一 个 特定 的 
Wi, SH ISB 为 零 ， 二 是 令 整 个 测 站 网 的 15B 和 为 零 。 

目前 ，iGMAS 上 海天 文 台 分 析 中 心 ， 采用 上 述 第 二 种 约束 方法 对 BDS/GPS, Galileo/GPS, 
GLONASS/GPS W ISB 进行 解 算 。 图 5、 图 6、 图 7 分别 表示 2018 年 9 月 17 日 到 2019 年 

9 月 19 日 10 个 相同 测 站 的 BDS, Galileo, GLONASS 相对 于 GPS 系统 伪 距 ISB 值 的 时 间 
序列 图 。 可 见 同一 时 间 段 内 同一 个 测 站 下 ，Galileo/GPS 和 GLONASS/GPS 之 间 1SB 非 
常 稳定 ，BDS/GPS 的 稳定 性 稍 差 。 
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5 BDS/GPS ISB 时 间 序 列 
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图 6 Galileo/GPS ISB 时 间 序 列 


有 关 多 模 GNSS 定位 的 研究 相对 较 多 ， 主 要 包括 相对 定位 与 单 点 定位 两 种 方式 ，ISB 
的 处 理 策略 也 略 有 差异 。 
在 相对 定位 中 ， 时 间 基 准 偏差 会 在 双 差 处 理 中 消除 ， 但 是 相对 的 接收 机 设备 时 延 还 会 存 
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7 GLONASS/GPS ISB 时 间 序 列 


在 ， 所 以 仍然 存在 系统 间 偏 差 ”。 在 相对 定位 中 对 ISB 的 处 理 又 分 为 两 种 情况 : (1) SABE 
间 频 率 相同 时 ， 采 用 紧 组 合 模 型 进行 解 算 ;使 用 该 方法 的 模糊 度 还 会 保持 整 周 特性 ”， 但 
是 需要 两 个 系统 间 有 相同 频率 的 信号 ， 相 当 于 将 两 个 系统 看 成 单 系统 进行 处 理 。(2) 当 系 统 
闻 的 频率 不 同时 ， 采 用 松 组 合 模型 ， 该 模型 分 别 对 不 同 的 GNSS 系统 独立 进行 解 算 ， 系 统 
间 不 进行 交叉 处 理 ， 使 用 该 方法 的 伪 距 和 相位 ISB 均 可 消除 ”。 由 于 GPS 和 Galileo 系统 
有 频率 相同 的 信号 (L 1/E1, L5/E5a)， 这 两 个 系统 常常 采用 紧 组 合 的 方式 进行 数据 融合 ” 。 
而 在 GPS/BDS 联合 处 理 中 ， 较 多 采用 松 组 合 的 相对 定位 模型 王国。 在 关于 GPS/Galileo 
联合 进行 相对 定位 的 TSB 研究 中 ， 发 现在 接收 机 相同 时 1SB 值 差异 很 小 ， 不 需要 考虑 ， 而 
接收 机 不 同时 伪 距 ISB 的 相差 很 大 ， 有 的 甚至 达到 几 百 纳 秒 ”。 

在 单 点 定位 中 ，ISB 可 作为 一 个 附加 参数 进行 计算 ， 但 是 未 知 参数 比较 多 ， 一 般 会 引 
入 一 个 约束 条 件 ”。 不 同 的 约束 条 件 ， 获 得 的 TSB 不 同 。 在 静态 单 点 定位 中 ， 一 般 使 用 精 
密 星 历 及 卫星 钟 差 处 理 ， 这 时 TSB 也 不 含 时 间 基 准 差 。 当 进行 动态 定位 时 ， 经 常会 采用 
广播 星 历 进行 数据 处 理 ， 算 出 的 1SB 包含 时 间 基 准 差 ”。 对 于 精密 单 点 定位 来 说 ， 相 位 
ISB 与 模糊 度 无 法 分 离 。 在 采用 模糊 度 浮 点 解 的 情况 下 ， 相 位 1SB 与 模糊 度 会 一 起 解 算 出 
kT, 

虽然 在 不 同 应 用 中 对 ISB 处 理 方式 有 差异 ， 但 很 多 研究 均 认 为 接收 机 类 型 和 天 线 类 
型 对 ISB 值 有 影响 ， 在 短期 内 ISB 值 稳定 性 较 好 一 。 因 此 ，Jiang™”、 张 辉 候 等 人 还 对 
ISB 进行 建 模 ， 根 据 一 周 的 ISB 值 预报 一 天 的 1SB 值 ， 并 进行 检验 对 比 ， 证 明 预 报 结果 
较 好 ， 具 有 一 定 的 可 预报 性 。 
目前 ，7SB 主要 还 是 作为 多 模 GNSS 数据 处 理 的 一 个 衍生 产品 ， 很 少 有 学 者 对 其 进行 
专门 的 研究 。 当 需要 进行 多 系统 联合 处 理 时 ， 一 般 根据 自己 对 系统 间 偏 差 的 理解 进行 处 理 ， 
还 未 真正 形成 一 个 完全 一 致 的 概念 。 而 且 当 采用 不 同 的 约束 条 件 ， 也 会 解 算出 不 同 的 ISB, 
限制 了 用 户 对 TSB 的 使 用 。 对 于 动态 定位 来 说 ， 如 果 能 够 将 7SB 进行 预报 ， 将 极 大 地 提高 
解 算 效 率 。 但 是 现在 对 于 ISB 的 研究 较 少 ， 也 极 少 有 人 对 长 期 的 1SB 进行 特性 分 析 。ISB 
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是 多 模 GNSS 联合 应 用 必须 考虑 的 一 个 偏差 ， 如 果 能 够 更 好 地 认识 15B， 有 利于 多 系统 数 


据 融 合 处 理 。 


6 总 结 


本 文 主要 介绍 了 四 种 在 GNSS 数据 处 理 中 与 设备 时 延 相关 的 1 


i 差 ( 码 间 偏差 、 未 校准 的 


设备 时 延 、 频 间 相 位 偏差 及 系统 间 偏 差 ) 的 定义 、 计 算 方法 和 研究 现状 等 。 讨 论 了 这 些 偏差 


的 来 源 及 其 对 GNSS 应 用 的 影响 ， 并 详细 介绍 了 在 实际 应 | 


] 中 每 个 偏差 的 具体 处 理 方法 。 


码 间 偏 差 可 认为 是 两 个 信号 伪 距 设备 时 延 组 合 。 当 使 用 两 个 信号 时 ， 伪 距 观测 量 必然 会 
受到 与 信号 和 频率 相关 的 DCB 的 影响 ， 因 此 在 观测 中 需要 考虑 它 。 
善 ， 在 许多 应 用 中 使 用 。 未 校准 的 相位 时 延 主 要 是 为 了 固定 整 周 模糊 度 而 提出 的 ， 一 般 研究 


的 是 划 


随 着 多 GNSS 的 迅速 发 展 及 应 月 


小 数 部 


目前 DCB 发 展 最 为 完 


分 FCB， 主 要 是 由 于 相位 设备 时 延 与 模糊 度 无 法 完全 分 离 。 为 了 更 快 更 好 地 
利用 三 频数 据 ， 以 其 中 一 个 卫星 钟 差 为 参考 ， 其 他 钟 差 再 利用 频 间 相位 偏差 与 参考 钟 差 进 行 
转换 。 系 统 间 偏差 则 是 在 多 个 系统 组 合 应 用 时 ， 由 于 不 同 GNSS 系统 基准 、 信 号 体制 等 差 
异 产生 的 系统 性 偏差 ， 是 与 设备 时 延 有 关 的 相对 偏差 。 


需求 的 不 断 提高 ，GNSS 相关 偏差 的 研究 越 来 越 受 到 


重视 ， 这 些 与 设备 时 延 有 关 偏 差 的 标定 、 监 测 和 合理 的 处 理 方法 ， 对 提高 GNSS 的 高 精度 


服务 性 


关 的 偏差 进行 系统 的 整合 归纳 ， 将 它们 其 体 地 定义 并 对 它们 
地 便利 今后 在 相关 邻 域 的 研究 。 
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Research Progress on GNSS Bias 


WANG Han-yu!?, SONG Shu-li!, ZHOU Wei-li!, CHEN Qing-ming? 


(1. Shanghai Astronomical Observatory, Chinese Academy of Science, Shanghai 200030, China; 2. 
University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: A review of the research conducted until present on the subject of Global Nav- 
igation Satellite System (GNSS) device time delay is here provided. In the GNSS data 
processing, due to the phase center of the antenna at the satellite and receiver is inconsisten- 
t with the signal generation point and reception point, device time delay will be generated in 
the course of the signal transmission process. So, device time delay may affect the accuracy 
of positions and other quantities derived from the observations. Generally, device time delay 
is highly relevant to the clock bias and ambiguity, which makes it difficult to completely 
separate it from them. Therefore, device time delay and its related factors can turn into 
the relative device time delay according to actual needs. If these relative delays can be bet- 
ter classified and distinguished, it will have important reference significance for the future 
high-precision GNSS applications. Several different deviations based on device time delay 
are introduced respectively: differential code bias, uncalibrated phase delays, inter-frequency 
clock bias and inter-system bias. The origin of deviations is discussed along with their effect 
on GNSS application, and descriptions of how delays can be estimated and in other ways 


handled in the data processing are provided. 


Key words: differential code bias; uncalibrated phase delays; inter-frequency clock bias; 


inter-system bias 


